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Abstract 

We study the macroscopic quantum tunneling , self-trapping phenomena in 
two weakly coupled Bose-Einstein condensates with periodically time- varying 
atomic scattering length. 

The resonances in the oscillations of the atomic populations are investi- 
gated. We consider oscillations in the cases of macroscopic quantum tunnel- 
ing and the self-trapping regimes. The existence of chaotic oscillations in the 
relative atomic population due to overlaps between nonlinear resonances is 
showed. We derive the whisker-type map for the problem and obtain the esti- 
mate for the critical amplitude of modulations leading to chaos. The diffusion 
coefficient for motion in the stochastic layer near separatrix is calculated. The 
analysis of the oscillations in the rapidly varying case shows the possibilty of 
stabilization of the unstable ir - mode regime. 
PACS numbers: 03.75 Fi, 05.30 Jp 
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I. INTRODUCTION 



Recently much attention has been paid to the investigation of macroscopic tunneling 
phenomena between coupled Bose-Einstein condensates (BEC) [OLE). The theory has been 
developed for weakly as well as for strongly overlapped condensates 0- ||. As a result, 
it was shown that, in oscillations of the relative atomic population, periodic processes of 
two types can exist. One which is characterized by a zero mean atomic inbalance function 
(< z{t) >= 0), corresponding to the periodic flux of atoms between two condensates due 
to the overlaps of the wavef unctions. The second type has a net atomic population balance 
< z(t) >t^ 0, corresponding to the so-called macroscopic self-trapping , with localization of 
atoms in one of the condensates, and a periodically varying population around a constant 
value. The different regimes are connected with the value of the initial phase differences 
9 = <pi — <p 2 between the condensates and the effective nonlinearity A ~ (ai + 0:2) N, where 
«j are the nonlinearity parameters depending on the atomic scattering length a s (which in 
turn depends on the two-atoms interactions), and N is the total number of atoms. 

If a periodic variation of the parameters of the trap and of the BEC is introduced, new 
phenomena show up. Examples of this are the resonant phenomena with a periodically 
varying trap potential considered recently for the case of weakly coupled BEC in || and for 
the strongly overlaped BEC in |TD[. The optical analog of this process is the electromagnetic 



wave transmission in a nonlinear coupler with periodic variation of tunnel coupling |TT],[T2 



Further interesting cases occur when a time-dependent atomic scattering length in BEC is 
considered. This is made possible due to the fact that a s depends resonantly on the external 
magnetic field or on the laser field |13|- ||18|| . This means that the scattering length can 
be varied in time, in particular, periodically. The main interest now is to investigate the 
influence of this periodic variation of atomic scattering length on the macroscopic interference 
phenomena. The central point here consists in that, as the interference phenomena are 
time-periodic processes in a nonlinear system, we can expect a nonlinear resonant response 
of the system under time-periodic perturbations, leading to a variety of phenomena, such as 
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hysteresis or chaos. 



II. FORMULATION OF PROBLEM. SLOWLY TIME- VARYING ATOMIC 

SCATTERING LENGTH 

The problem of weakly coupled BEC's in a double-well trap potential with time- 
dependent atomic scattering length a s can be described by the two-modes model 

ih^ = (E 2 + a 2 (t)tp 2 \ 2 )ij 2 - Kfa, (1) 

where the parameters E^a^K are defined by the overlaps integrals of the time dependent 
Gross-Pitaevsky eigenfunctions and derived in Pj. In particular, are the parameters 

characterizing the nonlinear interactions between atoms and a« ~ go = 47rh 2 a s (t)/m, where 
a s (t) is the time dependent atomic scattering length. This system of equations is valid in 
the approximation of a weak link between condensates. As comparison with the numerical 
solution of GPE shows that it is a good aproximation for z < 0.6. (see the definition of z(t) 
below). Anther limitation is connected with the time dependence of a s . For a harmonic 
modulation in time we should require that the perturbation does not introduce transitions 
between the ground state and excited states in the trap, whose energy gap is of the order 
of huo, where ujq is the harmonic frequency of trap. Then we have the restriction Q <C ujq. 
Another characteristic frequency is given by the coupling between condensates and is defined 
as u>l = 2K. For the weak coupling case when uo^ ujq the regimes of resonant Q ~ uol 
and rapid modulations of a s can be realized. The periodic modulation of nonlinearity in 
the Gross-Pitaevsky equation can also lead to the parametric instability of the collective 
excitations of the condensate [19|,20[|, that, in turn, leads to the break down of the two- 



modes approximation. Because these frequencies are far from the ones considered below we 
wiil use the model ([!]) in our analysis. 

Introducing new variables ipi = a/A^ exp(i6 , i ), z = (Ni — N 2 )/N T , N T — Ni + N 2 , $ = 
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#1 — 6*2, where iVj, ^ are the number of atoms and phases in the i-th trap, we get the system 



z t = — Vl - £ 2 sin$, (2) 
$ t = vA(t)z+ j= — - cos $ + AE , (3) 

V 1 — 2 

where A = (ai(t) + a2(t))N T /2K, t = uj^t, u L = 2K, u L being the frequency of the 
linear Rabi oscillations, AE = AE/2K, and v = ±1 for the positive and negative atomic 
scattering length, respectively. 

The Hamiltonian of the unperturbed system (a s (t) = const, AE = const, rj = 0) is 

vAz 2 

H = — cos($) + AE z. (4) 

Below, if not specified otherwise, we will consider the case of positive scattering length 
v = l,and AEq = 0. Consider the initial phase difference to be zero, $(0) = - the zero- 
phase mode. Depending on the number of atoms, N, and two-atoms interaction, a, i.e. A, 
there may exist a macroscopic quantum tunneling regime with < z(t) >= 0, corresponding 
to A < A c and a self-trapped regime with < z(t) >^ 0, corresponding to A > A c . The same 
is valid for the initial phase difference, $(0) = n - the 7r-phase mode. 

Let us consider now the case of periodic modulations of the atomic scattering length, 
when a s (t) = a (l + e cos fit), inducing a variation of A, i.e, A = A + Ai cos(Ot). Such a 
variation can be achieved by the variation of the external magnetic field [18| - a magnetic 



field induced Feshbach resonances for example. Near the resonance the scattering length is 
varied dispersively and can have positive and negative values. We will consider situations 
when we are not close to resonance and variations of scattering length are small, that is, 
Ai/Ao = e < 1. 

A. Resonances in zero phase modes 

Let us first consider the case where $(0) = 0. Take the case of z 2 1 and Ai > z. 
Taking into account that z t ~ $ ~ z 2 , assumed to be small, we obtain the equation for the 
relative atomic population 



z a + (A(f) + 1)2! - M±i^ = 0. (5) 
where, A(t) = A (l + esin(fit)). For e < 1 we can reduce the equation to 

z tt + - f3z 3 + eA cos(Qt)z = 0, (6) 

where uq = \J1 + A is the frequency of linear oscillations, and (3 = (A + l)/2. It is useful 



to introduce a new variable y — y (A + 1)2 < 1 and a scaled time r = tio . Then we have 
the equation 

Vrr + y- Piy 3 + e x cos(^ir)y = 0, (7) 

where (3 X = l/(2y/A + l),£i = eA /(A + l),^i = n/w . 

Let us consider the dynamics in the parametric resonance region when Qi = 2 + A. Using 
the results of perturbation theory ||21|| , we find that the parametric resonance occurs when 
the condition for the square of increment s 2 > is satisfied, where 

* 2 4[(y(l-^)) 2 -A 2 ]. (8) 

We have included here a nonlinear correction to the frequency 3/?ia 2 /8. Then we obtain the 
limit for the amplitude for parametric resonance to occur, 



8( ei - 2A) 
:W ; ., 



and z p = a p J vAo + 1. 

In Fig.(l) we present the oscillations of z(t) when the parameters are in the region of 
the parametric resonance, by integrating numerically Eq.(pj). 



B. Resonances in tt - phase modes 

We now come to the case 0(0) = tt, with for A < 1 and z 2 << 1. The system (|2|) is 
simplified and reduces to one equation for $, having the form of the equation of a double 
pendulum. Inspection of this equation shows that there exists a valley in the effective energy 
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V($>) around $ = 7r. The maximum of the valley depth is achieved when Ao — > 1. Thus 
we can search the solution of the system in this region of parameters assuming 2; C 1,$ = 
71 + S(t). It results that we obtain the following equation for z(t) 

z tt + (1 - A )z - Ai cos(nt)z = 0. (10) 

Here Ao < 1 .The parametric resonance in the oscillations of the atomic population appears 
when Q = 2vT^A~o + A. The parametric instability occurs when Ai/ (2a/T^-~Ao) > A. The 
increment s is, in this case, 



For example, when A = 0.36, A = 0.1 we have Q = 1.7. In Fig. (2) we plot the oscillations 
of z(t) for this choice of parameters, and with Ax = 0.162, by a nummerical integration of 
Eq.(§). 

Note that the frequency of linear Rabi oscillations is ujl = 2K, or one in the units used. 
Thus the resonance in n -mode occurs at frequencies lower than the Rabi frequency. It is 
worth mentioning that the metastable 7r-phase mode has been recently observed in the weak 
link separating two reservoirs of superfluid He 3 W\ . 



III. CHAOTIC DYNAMICS IN OSCILLATIONS OF THE RELATIVE ATOMIC 

POPULATION 

The oscillations induced by the periodic perturbations of the scattering length may be 
comples and become chaotic, for certain regions of the parameters. Note that the unper- 
turbed system is equivalent to the Duffing oscillator. Thus, as is known, the resonance 
overlaps under periodic perturbations can lead to chaotic variations of z(t) p3|-p5|. In such 



cases, to study the chaotic motion, it is useful to look for the motion near the separatrix of 
the unperturbed system. The total Hamiltonian is H = Hq + V, where Hq is unperturbed 
Hamiltonian (4) and the perturbation is V = Ai cos(f2t)z 2 /2. The value of the Hamiltonian 
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on separatrix is H s = 1. The separatrix divides the regions with macroscopic quantum 
tunnelling, with < z >= 0, from the self-trapped regions where < z >^ 0. The expressions 
for the separatrix solution are: 

z 8 (t) = iJ^sech(V2at) , 



a 2 sech. 2 (y/2at) tanh 2 (v / 2at) 
2K 2 (b - asech^v^t)) 



sin = nT ^ t L -—,2/ 7^ . ( 12 ) 



where a = 1 1 — A|/2, b = A 2 /8. For example, for the zero-phase mode the separatrix is given 
by the hyperbolic fixed points $ = ±7r, z — 0. 

The existence of chaos can be proved by calculating the Melnikov function. The calcula- 
tions are similar to the ones performed in and shows that the Melnikov function has an the 
infinite number of zeros, so the existence of chaotic regimes in the atomic population oscilla- 
tions is proved. In Fig. (3) we plot the Melnikov distance versus the frequency Q for A = 10. 
We see that the maximum value is achieved for Q « 3.89. In Fig. 4 we plot the difference of 
atomic population z(t) as function of time for A = 10, z(Q) = 0.4, Ai = 0.2, Q = 3.89. 

It is interesting to obtain an analytical estimate for the critical amplitude of modulation 
Ai leading to the chaotic behavior in the relative atomic population. Using the expression for 
the separatrix (|12"D we can calculate the Melnikov- Arnold integral and find that the energy 
change for a perturbation Ai cos(f2t + ip) (where ip is the phase) is 

AH - 1 Jit + v » dt a - aJ^W ' (13) 

where [...] denotes the Poisson bracket. 

The period of the unperturbed motion is T = AkK(k)/(CA), where ||: 

n 2 = -[l+ , H ° A ~ 1 ] , (14) 
2 [ VA 2 + 1 - 2H A l ' 1 ; 



C 2 = ^ [(# A - 1) + VA 2 + 1 - 2# A] , (15) 

where K(k) is the complete elliptic integral of the first kind. Near the separatrix, when 
H — > 1, we have k 2 — > 1. Introducing the parameter 5 = 1 — H Q , H < 1, 5 << 1, and 



taking into account that K(k) — > ln(4/Vl — ^ 2 ), when k 2 — > 1, we obtain the estimate for 
the period of oscillations near the separatrix 



2 /4*/2(A-r. 
, In ( V ^_ 1. (16) 



Using the expressions fllTf ) and ([II]) we find the whisker-type map for our problem 



S n +i = $n + asin(^. 



n/i 



Tpn+1 = Ipn + 7 l n f 7t= I • ( 17 ) 



+ 1, 



where 7 = 2Q/\/A^T and // = 4^/2(A - 1)/VA. 

This map can be simplified, using the linearization around the fixed points [21|. The 
fixed points are defined by 

27r/ = 7 ln(^=] , / = 1,2... (18) 



Let us introduce the dimensionless energy I n , defined by 

I n = -^(5 n -5f). (19) 
°f 

Substituting this expression into Eq. (|TTD , and redefining I n — > I n /2, we obtain the stan- 
dard map 

4+i = In - Ksm(ip n ) , 

tpn+1 = VVi + In+1 , (20) 

where the parameter K is 

K=^. (21) 
25f 

As is well known, chaos appears when K > 1 (exact value is i^* = 0.9716). Then we have 
the estimate for the critical amplitude of modulation 

> 4° A 3 sinh(7r^/2 v / A^T) 
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From the above, we may calculate the diffusion coefficient for the motion in the stochastic 
layer 

D = <( A f>. (23) 

It results that we have D = K 2 /2. For frequencies Q < 2\/A — 1/tt we have the estimate 
„ 2 1( W(A-l) 3 ft 4 , 47t/ v / A^T, 

D ~ ^ ' exp( ^ ) . (24) 

For £ = 0.1, Q = 2.5, A = 9, I = 1, the time of exit from the stochastic layer is 10 i.e. ~ 4 
periods of the oscillations of the atomic scattering length. 



IV. RAPIDLY VARYING SCATTERING LENGTH 

In the case when the frequency of oscillations of the atomic scattering length is much 
larger than the tunneling frequency ujl = 2K. i.e. A = A(t/e), where e = ujl/Q « 1, 
it is useful to derive an averaged set of equations. Using a multiscale approach |26[] and 
expanding z = z + ez\ + $ = $ + e$i.... we can derive the system of equations for the 
slowly varying functions z, The limits of validity is z 2 << 1. 

The averaged system turns out to be: 



Zt = -Vl - £ 2 sin($)(l - Sz 2 ), (25) 
§ t = u<k>z+ * cos ^) (i + 25 - 35z 2 ). (26) 

where 5 is proporttional to e 2 , and for the harmonic modulation is Af/(4f2 2 ). The corre- 
sponding Hamiltonian is 

v < A > 



H = z 2 - v / T 3 i 2 cos(<l)(l - 5z 2 ) . (27) 

With respect to the unperturbed case, a new stable point can occur. A bifurcation is 
observed when 5 is increased. The unstable point (in the undisturbed case) at z = 0, $ = n 
may become stable under rapid perturbations. This phenomenon is the analog of the Kapitza 
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stabilization of the unstable fixed point of the pendulum by rapidly oscillating the suspension 
point. By passing though the bifurcation, a change in the topology of the phase portrait 
happens. The value of S for this to occur depends on < A >. To have an estimate of the 
critical value of 5, consider z 2 « 1 and $ = tt + ip, ip << n. The equation for z reads: 

z tt - (< A > -1 -2S)z = . (28) 

Thus, when 5 > 5 C = (< A > — 1)/2 we have oscillatory motion around (0, n) in the phase 
plane. When 5 < 8 C we have an hyperbolic point, and (0, tt) is no longer stable. For 
example, with A = 2, the critical 5 is 0.50. For < A >= 1, the point is always stable, 
as S c = 0. To support these results, we have numerically integrated the original systems 
given by Eq.(|2]) with an explicit harmonic fast-time variation for A(t). The above described 
features immediatly show up. 

Note also that we can cross the sign of a s from positive to negative. Experiments with 
very a large monotonic change of sign show that the BEC shrinks very fast. It would be 
interesting to check if the condensate with attractive interaction of atoms stable under the 
rapid variation of scattering length. This problem requires, however, separate investigations. 
In any case in quasi ID geometries collapse is suppressed and this analysis is relevant. 

In order to better understand the physical consequences of the averaging, we note that 
for small z 2 « 1 we can rewrite the Hamiltonian in the form 

^=^-/^co S ($), 

where z T = (1 + 5)z, A r = (1 — 25) < A >. Thus we can conclude that the result of the 
averaging consists, for a fixed initial phase difference between condensates, in the lowering 
of the critical relative population for self-trapping regime to occured and in the increase of 
the critical value of the nonlinear parameter A. We come to the picture of a more rigid 
pendulum in comparison with the constant case. The threshold for the switching from the 
macroscopic quantum tunneling regime to the self-trapping regime is shifted to lower initial 
values of the atomic imbalance. 
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V. CONCLUSION 



In this paper we have studied the new effects occuring due to time-periodic variations of 
the atomic scattering length. The resonances in oscillations of the atomic imbalance for the 
different regimes have been analyzed. 

The interaction of resonances occuring under oscillating scattering length can lead to 
complicated behavior of the atomic oscillations and gives rise to the phenomenon of chaotic 
macroscopic quantum tunneling. In this work we prove the existence of such regime, using 
the Melnikov function approach and calculating the Melnikov distance, characterizing the 
width of the stochastic layer near the separatrix of the unperturbed system. We derive a 
whisker type map for the problem and obtain the estimate for the critical value of amplitude 
A x leading to the chaotic motion. The diffusion coefficient for the motion in the stochastic 
layer is calculated and shown that the time to cross the stochastic layer is the order of 
few periods of modulations. We also consider the evolution of the system under rapidly 
varying oscillations of the scattering length and derive the averaged system for the slow- 
time variations of the relative atomic population z(t) and the phase $(£). The analysis of 
the fixed points shows that the stabilization of the system under rapid perturbation, when 
the unpertubed is unstable, is possible in the n -phase mode. 
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FIG. 4. The chaotic oscillations of the relative atomic population z(t) for 
n = 3.89, A = 2, A = 0.2, $(0) = 0, z(0) = 0.6 
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